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The possible functional forms of the effective conductivity cr^f f of the randomly 
inhomogeneous two-phase system at arbitrary values of concentrations are dis- 
cussed. A new functional equation, generalizing the duality relation, is deduced 
for systems with a finite maximal characteristical scale of the inhomogeneities 
and its solution is found. A hierarchical method of the construction of the model 
random inhomogeneous medium is proposed and one such simple model is con- 
structed. Its eff'ective conductivity at arbitrary phase concentrations is found in 
mean field like approximation. The derived formulas (1) satisfy all necessary in- 
equalities and symmetries, including a dual symmetry; (2) reproduce the known 
formulas for 0-^/ / in weakly inhomogeneous case. It means that in general 
Ue// of the two-phase randomly inhomogeneous system may be a nonuniver- 
sal function and can depend on some details of the structure of the randomly 
inhomogeneous regions. The percolation limit of the randomly inhomogeneous 
two-phase systems is briefly discussed. 



The electrical transport properties of the disordered systems have an impor- 
tant practical interest. For this reason they are intensively studied theoretically 
as well as experimentally. In this region there is one classical problem about 
the effective conductivity agf / of inhomogeneous (randomly or regularly) het- 
erophase system which is a mixture of N{N > 2) different phases with different 
conductivities ai,i — 1,2, N. We confine ourselves here by the simplest case 
of the two-dimensional heterophase systems with N = 2. Despite of its rela- 
tive simplicity only some few general exact results have been obtained so far. 
Firstly, there is a general expression for cTg/ / in case of weakly inhomogeneous 
isotropic medium, when the conductivity fluctuations Sa are smaller than an 
average conductivity (a) {Sa (a) ) similar to the analogous expression for 
the permitivity of the corresponding dielectric mixture 

--<">('-^)' « 

where D is a dimension of the system. In our simplest case of two-dimensional 
two-phase system {a) = xai + (1 — x)a2, (cr^) — (c)^ = x{l — x){Aa)'^, where 
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a; is a concentration of the first phase, Act = cti — CT2 , and the first formula 
takes a form 



,Act xfl - .t)(Act)2 



/ , , X Act x(1 - x)(Act)2 , ^ ,^ 

ct+(i + (:c-1/2)- ^ 4 



' + 

where ct+ = (cti + CT2)/2. 

Another general formula has been obtained for the case of a small concen- 
tration of one phase (for example, with a conductivity CT2 ) [Q 

CTe,,=CT,(l-(l-x)fc^), (2) 
\ CTi + CT2 / 

where 1 — a; ^ 1 is a small concentration of the second phase and a round form 
of the inclusions of this phase is suggested. 

The further progress in the solution of this problem is connected with papers 
[|[ ^ . It was shown there that there is a dual transformation transforming one 
phase into another. This transformation allows to find an exact formula for 
(Jeff in case of systems with equal concentrations of the phases a; = a;c = 1/2 

i 

(Jeff = ^/(Jl(y2■ (3) 

This remarkable formula is very simple and universal since it does not depend 
on the type of the inhomogeneous structure of the two-phase system. For sys- 
tems with unequal phase concentrations the dual transformation gives a relation 
between effective conductivities at adjoint concentrations x and 1 — x or in 
terms of a new variable e — x — Xc ( ^1/2 < e < 1/2 ) at e and — e 

CTe//(x, CTi,CT2)CTe//(l - a;,CTi,CT2) = CT1CT2 = CTg// (e, CTi , CT2 )CTe// CTi , CT2 ) . (4) 

Eq.(4) means that a product of the effective conductivities at adjoint concen- 
trations is an invariant. Due to this relation one can consider (Jeff only in the 
regions a;>a;c (e>0)ora:<a;c (e<0). 

In the following papers this dual transformation was generalized on systems 
in magnetic field ||, ^, ^ , polycrystals ||^ and heterophase (with > 2 ) 
systems |,|. The formula (4) was also checked and improved for some regularly 
inhomogeneous systems (of the chess-board type) jl^, |l^ . 

However the main interest in this problem has a formula for the effective con- 
ductivity at arbitrary phase concentrations. Another question appears naturally 
in this case: is this formula also universal or it can depend on the structure of the 
two-phase system? We will give here a general analysis of the possible functional 
forms of the effective conductivity in the case of arbitrary phase concentrations 
and will show that, due to the dual symmetry, in some cases the corresponding 
effective conductivities can be found in the explicit form. Then we will intro- 
duce the finite maximal scale averaging approximation (FMSA-approximation), 
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deduce in the framework of this approximation a new equation, connecting the 
effective conductivities at different concentrations, and find its solution and its 
physical meaning. We will introduce also a new model of the random two-phase 
system which gives in the same approximation an effective conductivity, coin- 
ciding with one obtained by series expansion. All these results demonstrate 
that, in general, a formula for the effective conductivity may be nonuniversal 
even in the two-phase case. At the end of the paper we will briefly discuss some 
peculiarities of the percolation limit. 

Let us start our investigation of the isotropic classical random two-phase 
system in the case of arbitrary concentrations with a general analysis. The 
basic equation is the Ohm law, connecting the local current j(r) and local 
electric field E(r), 

j(r) = a(r)E(r), (5) 

where a(r) is a local conductivity, ft must be supplemented with the corre- 
sponding boundary conditions on the boundaries of two phases ||T[ 

jI^jI El^El (6) 

where n,i denote the normal and tangent components and 1,2 correspond to 
different phases. The electric field E(r) is a curl-free field 

V X E(r) = 0, (7) 

and the current field is a divergenceless field 

V-j(r)=0. (8) 

The effective conductivity CTe// of the isotropic system can be defined as a 
proportionality coefficient between averaged values of j and E over the area 
of the system 

j = f^eZ/E, 

]=Jj{v)dr/S, E = jE{v)dr/S (9) 

where S is an area of the system. Due to the linearity of the defining equations 
an effective conductivity of the random systems must be a homogeneous function 
of degree one of ai,i = 1, N. In our case N — 2 and it is convenient to use 
instead of ai,i = 1,2 another combinations of ai : a± = (fxi ± cr2)/2. Then, 
introducing a new variable z = j a^^ (—1 < < 1), an effective conductivity 
can be represented in the following, symmetrical relatively to both phases, form 

cre//(e,cr+,(T_) = cr+/(e, (T_ /cr+) = (T+/(e,z), (10) 

where f7e//(e, f^, f7_) and f{e,z) must have the next boundary values 

t^e//(l/2,Cr+,(T_) = (Ti, cre//(-l/2,cr+,(T_) = (72, 

/(l/2,z) = H-z, /(-l/2,z) = l-z, /(e,0) = l. (10') 
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The duality relation, being a consequence of a duality of gradient and tangent 
fields in two dimensions js), takes in these variables the form 

f{e,z)f{-e,z)^l~z\ (11) 

from which it follows that at critical concentration 6 = 

/(0,z)-yT^. (3') 

Strictly speaking, this form of a duality relation is also a consequence of another 
exact relation for the effective conductivity, taking place at arbitrary concen- 
trations for systems with the similar random structures of both phases of the 
system, 

CTe//(e,Cri,cr2) = (7^// (-£, (72 , ^1 ) . (12) 

It means that the effective conductivity of the random two-phase system must 
be invariant under substitution of these phases ( (Ti < — > <T2 ) with the corre- 
sponding change of their concentrations x < — > \ — x (or e — > — e ). In the new 
variables it means that 

/(e,z)-/(~e,-z), f{-e,z)^f{e,-z). (13) 

For this reason a duality relation can be written also in the form 

f{e,z)f{e,~z) = l^z\ 

It follows from (13) that the even ( /s ) and odd {fa) parts of /(e, z) relatively 
to e coincide with the even ( ) and odd ( ) parts of /(e, z) relatively to 
z 

r (6, z) ^ 1 (/(e, z) - fie, -z)) = /,(6, z) ^ i (/(e, z) - /(-e, z)) , 

/^(e, z)^\ (/(e, z) + /(6, -z)) ^ fs{e, z)^\ (/(e, z) + f{-e, z)) . (14) 

The simplest way to satisfy (13) corresponds to the following functional form of 
f{e,z) 

f{e,z) = fiez,e^,z') (15) 

One can conclude from (15) that in this case an expansion of f{e.,z) near 
the point e — z = does not contain terms linear in e and z separately. 
Moreover, one can reduce to the form (15) any function satisfying (13), since 
every combination of functions of two variables, odd relatively to the inversion 
of their arguments, can be made even by multiplying or dividing them by ez 

F{-e, z) = -F(e, z) = F{e, -z), F(e, z) ^ ez$(e^ z^), 

$(e2,z2) = (ez)-iF(e,z). 
Analogously, the odd part fa can be represented in the form 

/,(e,z) = 2ez$(e,z) (16) 
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where $ is an even function of e and z (the coefficient 2 in front of ez is 
chosen for further convenience). 

The duahty relation (11) in general case is not enough for the complete 
determination of / , it only connects /q and /s 

-f! = l~ z\ (17) 

It means that /a and considered at fixed z as the functions of e satisfy to 
hyperbolic relation with a constant depending on z. The relation (17) allows 
to express /(e, z) through its even or odd parts 



/(6,Z) = + + 1 - Z2 ^ /. ± V/l - 1 + (18) 

For this reason it is enough to know only one of these two parts. Usually one 
prefers to choose an antisymmetric part as more simple. It follows from (1') 
that in the weakly inhomogeneous case the odd part has the simplest form (16) 
with <I> = 1 

/a(e,z) = 2ez. (19) 

As is well known, the substitution of (19) into (18) gives an effective conductivity 
in the effective medium (EM) approximation jj] 



cr^^^(e,z) = (7+ 2ez+ v/(2ez)2 + 1 - z2 , (20) 

This expression, being continued on arbitrary concentrations e = a; — 1/2 and 
inhomogeneities z, reproduces in the corresponding limits both formulas (1') 
and (2). 

Though at first sight a duality relation is not enough for finding / , usually 
systems with such symmetry have some additional hidden properties, permitting 
to obtain more information about function under question. Moreover, in some 
cases they can help to solve problem exactly (see, for example, Q). Having 
this in mind, we will try to investigate the duality relation in more details. Since 
a homogeneous limit z ^ is a regular point of / it will be very useful to 
consider a series expansion of / in z 

oo 

/(e,z) = ^/fc(e)zVfc!, (21) 



where due to boundary conditions (10') 

/o = 1, /i(e) = 2e. (22) 

For every fixed z (0 < z < 1) a function / must be a monotonous function of 
e. Substituting the expansion (21) into (11) one obtains the following results: 

(1) in the second order in z it reproduces the universal formula (1'), thus 
the latter can be considered as a consequence of the duality relation; 

(2) in general case there are the recurrent relations between /2fc and /2fc-i, 
corresponding to the connection (17); 
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(3) /2/c+i(e) are odd polynomials of e of degree 2k + 1 and /2fe(e) are even 
polynomials of e of degree 2k in agreement with (15). 

Taking into accoimt boundary conditions (10') and an exact value (3'), one 
can show that the cocfBcionts fk must have the next form 

/2fe+i(e) = e(l - 4e2)g2fe_2(e), k>l, 

/2fe(e) = (l-4e2)/i2fe-2(e), fc > 1, (23) 

where g2k-2 and h2k-2 are some even polynomials of the corresponding degree 
and free terms of /i2fe-2 coincide with the coefficients in the expansion of (3') 

\/l-^2 = 1 - z^/2 - z*/8 - zVl6 - zVl28 - 2^7256 + ... (24) 

It follows from (23) that /a is completely determined up to overall factor num- 
ber go. Since /4 is determined through lower fk (fc = 1,2,3) 

fi = ^hh - 3/1 = (1 - 4e2)[(85o + 12)e' - 3], (25) 

it is also determined by the coefficient qq. The case go = = g2k+i {k > 1) 
corresponds to the EM approximation and in this case /2fe(e) ~ (1 — 46^)*^ . 
Thus we see that in general case the arbitrariness of / is strongly reduced 
by boundary conditions and known exact value and that the third and fourth 
orders are determined only up to one constant. Tlicni one can suppose that 
any additional information about function / can determine this constant or 
even the whole function. For this reason one can ask: what kind of functions 
can satisfy the duality relation (11)? In order to answer on this question it is 
convenient in the case z 7^ 1 to pass from / to / = // -\/l — z"^. Then 

/(e, z)f{-e, z) = l = /(e, ^)/(e, -z). (11') 

The duality relation gives some constraints on the possible functional form of 
f{e,z). For example, assuming a functional form (15), one can write out a 
simple function: 

/(e, z) = exp{ez(j){e, z)), (26) 

where c6(e, z) is some even function of its arguments. It is easy to see that it 
automatically satisfies eq.(ir). Let us now consider two simple cases, when (a) 
4>{e,z) depends only on z, (b) it depends only on combination ez. Expanding 
the corresponding functions / in series one can check after some algebra that 
one can now determine all polynomial coefficients unambiguously! Another way 
to see this is to apply boundary conditions directly to the function (26). In the 
case (a) one obtains 

<^(z) = l/zlnl±^, /(e,z)=(^i±iy. (27) 

In case (b), when ^ depends only on combination ez , one finds 

1 l + 2ez - /^i±2ffV^' 
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Series expansions of (27) and (28) coincide exactly with the corresponding ex- 
pansions mentioned above. Unfortunately, analyzing an expansion of (26) under 
an assumption of regular behaviour of the coefficients ^fe(e) from a series ex- 
pansion of (j){e, z) 

oo 

0(e,z) = ^<^fc(e)02fc/^, 



in e 

oo 

<^fe(e) = 5^<^«(ef A!, 



one can show that now again the boundary conditions (10') (or (23)) cannot 
define all coefiicients completely. For example, /s and /4 contain one free 
parameter c6io : .9o =6((/)io — 1). It is very important to find any other solutions. 
Below we will present two simple models having the effective conductivity just 
of two forms found above. 

Firstly, we will obtain an additional equation for effective conductivity using 
a new composite method, which we will call as a finite maximal scale averaging 
approximation. In framework of this method one divides an averaging procedure 
on two steps. This approximation can be considered as some modification of 
the mean field approximation and it is applicable to the inhomogeneous systems 
with a finite maximal scale of the inhomogeneities. Below in this section we will 
use temporally, for clarity, the concentration variable x instead of e. Let us 
consider two-dimensional two-phase randomly inhomogeneous system. It is easy 
to see that, in general, its effective conductivity will depend on a scale I of a 
region over which an averaging is done. This takes place due to the possible 
existence of different characteristic scales in the inhomogeneous medium. In the 
most general case there will be a whole spectrum of these characteristic scales. 
This spectrum can be very different from discrete finite till continuous infinite 
and will define the inhomogeneity structure of the system. For this reason 
this spectrum can depend on concentration x. It is obvious that the effective 
conductivity of the system can depend on this spectrum as a whole as well as 
on which area an averaging is fulfilled over. Suppose, for the simplicity, that the 
randomly inhomogeneous structure of our system has the scale spectrum with 
a maximal scale lm{x), which is finite for all x in the region 1 > a; > 1/2 (or 
1 — a; in the region 0<1 — a;<l/2). Let us assume that we know an exact 
formula for aeff{x,z) of this system, which is applicable from scales I > Im- 
It means that this formula for apfj{x,z) takes place after the averaging over 
regions with a mean size I > Im and does not change for all larger scales I ^ Im- 
Now consider a square lattice with the squares of length II ^ Im and suppose 
that they have the effective conductivities corresponding to different values of 
the concentrations xi and X2 with equal probabilities p=l/2 (see Fig. 1) 

After the averaging over the scales I > II one must obtain on much larger 
scales I ^ II the same effective conductivity, but corresponding to another 
concentration x = {xi + X2)/2. This is possible due to the similar random 
structure of different squares and due to the conjectures that in this model: 
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(b) 



Figure 1: (a) An elementary square of the model with round inclusions, (b) a lattice 
of the model, the numbers 1,2 denotes squares with the corresponding concentrations. 



(1) there are only two maximal characteristic scales lm{xi) = h {i = 1,2) 
and the lattice square size II ^ li,lm{x) , (2) the averaging procedures over 
these scales do not correlate (or weakly correlate) between themselves. Thus 
for compatibility all concentrations must be out of small region around critical 
concentration Xc, where U or lm{x) can be very large. We will call further this 
set of the conjectures the finite maximal scale averaging approximation (FMSA 
approximation). It can be implemented for systems with compact inhomoge- 
neous inclusions with finite Im ■ From the other side the effective conductivity 
on scales I ^ II must be determined by the universal Keller - Dykhne formula 
(3). Thus we obtain the next functional equation for the effective conductivity, 
connecting aeff{x,z) at different concentrations, 

aeff{x,z) = ^Jaeff{xi,z)aeff{X2,z), X={xi+X2)/2, {Xi,Xj^Xc). (29) 

It must be supplemented by the boundary conditions (10'). This equation can 

be considered as a generalization of the duality relation (4), the latter being a 
particular case of (29) at xi + X2 = 1. It follows from (29) that, due to the 
exactness of the duality relation, it really works for all concentrations x, except 
maybe of small region near x = Xc and z = 1 (this region corresponds to the 
singular region of the percolation problem, see also below a discussion of the 
percolation limit). One can easily check that CTe// in low concentration limit 

(2) satisfies (29). Moreover, one can find an exact solution of this equation. It 
has an exponential form with a linear function of x (or e ) 

aef f{x, z) = aic:xjp{ax + b), (30) 

where the constants o, b can be determined from the boundary conditions 

a = —b, exp& = (T2/o'i. (30') 

Substituting these coefficients into (30) one obtains 



(Jeff{x,z) = ai (a-2/CTi) 
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The solution (31) satisfies all symmetry relations (4, 11', 13) and can be repre- 
sented in the exponential form 



(Te/ f{x,z) = cr+\/l — exp In J. (31') 

which exactly coincides with the case (a) from (27) with (^(e, z) — ^\vl . Its 
odd part has a form 

fa{e, z) = sinh(ez0(e, z)), (32) 

and satisfies (17). 

It is interesting to note that the form of the solution (31) means that in the 
considered approximation one has effectively an averaging of the log a since it 
can be represented as 

logCTe// = (log Cr) = a;l0gcri + (1 - .T)logC72. (33) 

This coincides with a note made in ||] for the case of equal concentrations 
a; = 1/2 and with analogous results obtained later for random system in the 
theory of the localization jl^. One can check that (31) reproduces in the weakly 
inhomogeneous limit the universal Landau - Lifshitz expression (1). In the low 
concentration limit of the second phase it gives 

I — z 

tTe//(a;,z) =(Ti(l + (1-x) log -— + ...), 1-2;<1, (34) 

J_ ~\~ z 

what coincides with (2) in the weakly inhomogeneous case. Note that the expan- 
sion (34) contains the coefficients diverging in the limit \z\ 1. Such behaviour 
of the coefficients denote the existence of a singularity in this limit (see below 
a discussion of this percolation limit). 

Now we will show that an existence of the different functional forms for the 
effective conductivity denotes a nonuniversal character of this value for two- 
phase systems with different random structures. We will construct a model 
of two-dimensional isotropic randomly inhomogeneous two-phase system, using 
the composite method introduced above, and find a mean field like expression 
for its effective conductivity Uef f{x,z) in case of arbitrary phase concentrations 
X. The derived formula satisfies again to all necessary symmetries, including a 
dual one, and coincides with the example case (b) from (28), realizing the second 
variant, when f{x, z) depends (in this approximation) only on one combination 
of variables ez and this dependence is described by the function analytical at 
small values of this variable. 

Let us consider the following two-dimensional model. There is a simple 
square lattice with the squares consisting of a random layered mixture of two 
conducting phase with constant conductivities ai,i = 1,2 and the correspond- 
ing concentrations x and 1 — a;. A schematic picture of such square is given in 
Fig.2. 

The layered structure of the squares means that the squares have some pre- 
ferred direction, for example along the layers. Let us suppose that the directions 
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(a) 



(b) 



Figure 2: (a) An elementary square of the model with a vertical orientation, the 
dotted regions denote layers of the second phase; (b) a lattice of the model, the small 
lines on the squares denote their orientations. 

of different squares are randomly oriented (parallelly or perpendicularly) rela- 
tively to the external electric field, which is directed along x axis. In order 
for system to be isotropic the probabilities of the parallel and perpendicular 
orientations of squares must be equal or (what is the same) the concentrations 
of the squares with different orientations must be equal = p± = 1/2. 

Such lattice can model a random system consisting from mixed phase re- 
gions, which can be roughly represented on the small macroscopic scales as 
randomly distributed plots with the effective "parallel" and "serial" connec- 
tions of the layered two-phase mixture (Fig. 2). The lines on the squares denote 
their orientations. 

This structure can appear, for example, on the intermediate scales when a 
random medium is formed as a result of the stirring of the two-phase mixture. 
The corresponding averaged parallel and perpendicular conductivities of squares 
(T| I (x) and cr_L {x) are defined by the following formulas 

fT||(a;) = X(7i + (1 — x)a2 = cr+{l + 2ez), 

Thus we have obtained the hierarchical representation of random medium (in 
this case a two-level one). On the first level it consists from some regions (two 
different squares) of the random mixture of the two layered conducting phases 
with different conductivities ai and a2 and arbitrary concentration. On the 
second level this medium is represented as a random parquet constructed from 
two such squares with different conductivities (T|| and (j± , depending nontriv- 
ially on concentration of the initial conducting phases, and randomly distributed 
with the same probabilities pi = 1/2 (Fig. 2). This representation allows us to 
divide the averaging process into two steps (firstly averaging over each square 
and then averaging over the lattice of squares) and implement on the second 
step the exact formula (3). This can be considered as some modification of 
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the FMSA approximation. As a result one obtains for the effective conductiv- 
ity of the introduced random two-phase model the following formula, which is 
applicable for arbitrary concentration 

(36) 



^eff{e-,z) = CT+Vl - z'^f{e,z), f{e,z) 



1 - 2ez 



This function has all necessary properties and satisfies equation (11') and sym- 
metry (13). It coincides with the type (b) from (28) and has another possible 
functional form, automatically satisfying the duality relation (11') 

f{e,z) = B{e,z)/B{-e,z) (37) 

with a function B{e,z) = [1 -|- 2ez]i , which depends only on the combination 
ez. 

It is interesting to compare this formula with the known general formulas. In 
order to do this one needs to find an asymptotic behaviour of the derived effective 
conductivity in different limiting cases. Let us consider firstly its behaviour for 
small phase concentrations. 

(a) In case of small concentration of the first phase x one gets 



(2xz \ 



(38) 



It follows from (38) that an addition of small part of the first higher conducting 
phase increases an effective conductivity of the system as it should be. 

(b) In the opposite case of small concentration of the second phase 1 — a; <C 1 
one obtains 

ae//(a=,.)^ai(l-^^^), (39) 

i.e. an addition of the phase with smaller conductivity decreases asff. It is 

worth to note that both these expressions for arbitrary values of the conductivi- 
ties cJi and (72 differ from equation (2) and coincide with it only in the weakly 
inhomogeneous case z <^ 1 . It must be not surprising because a form of the 
inclusions of the second phase in this model has completely different, layered, 
structure. In the low concentration expansion as well as directly in formula (36) 
one can see again that the divergencies appear in the limit 2; — > 1. 

(c) In case of almost equal phase concentrations a; = l/2-|-e, e<Cl one 
obtains 

/(e,0)~ l + 2e0, aeff{e,z):^a+VT^{l + 2ez). (40) 

The Keller - Dykhne formula (3) is reproduced for equal concentrations. 
The antisymmetric part of / has the following form 

~ Pie,z)-1 2ez 

/a(e,2) = ^r7^ ^ = 7; . 211/2 - (^1) 

2f{e,z) [1-Ae^z^]' 
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It follows from formula (41) that for e ^ 1 and (or) z ^ 1 the odd part fa 
coincides in the first order with the corresponding expression from the effective 
medium theory. The corresponding function $ is 



$(e,z) = ^^^^^. (42) 

In other words it is an analytic function of ez near ez = . Basing on this 
formula one can conjecture that an exact expression for fa{e,z) of this model 
will have a similar structure in general case with ai or 1 — z^^l. One 
must note that at the same time the formula (36) does not satisfy the equation 
(29) except of the trivial case xi = X2- 

For the comparison of the different expressions of the effective conductivity 
we have constructed three-dimensional plots of /(e, z) in the EM approxima- 
tion, in the FMSA approximation and of the hierarchical model in FMSA-like 
approximation (fig. 3, 4). 

It follows from these plots that all three formulas for a^f / , despite of their 
different functional forms, differ from each other weakly for z < 0, 8 due to 
very restrictive boundary conditions (10') and the exact Keller-Dykhne value. 
This range of z corresponds approximately to the ratio o^jox ~ 10~^. For the 
smaller ratios a difference between these functions become distinguishable. 

Now let us consider in the more details the derived formulas for aeff{e, z) 
in case when a2 0{z 1). It is clear that for regularly inhomogeneous 
medium one can always construct such distribution of the conducting phase 
that f7e//(e,l) will differ from zero for all 1/2 > e > —1/2. But in the case 
of randomly inhomogeneous medium the limit (T2 — > is equivalent to the well 
known percolation problem [14,15]. In terms of z it corresponds to the limit 
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z 1 and is also similar to the superconducting limit ai ^ oo . Strictly 
speaking, an implementation of the duality transformation (4) is not obvious in 
this case. However, if one supposes that the dual symmetry relation (4) fulfills 
in this limit too due to a continuity then it follows from (4) that 

^effie)aeff{-e) = 0. (43) 

The relation (43) does not contradict to the known basic results of the per- 
colation theory that Ce// (e) = for e < and Ce/ / (e) ^ for e > . 
Moreover it follows from the general formula (17) that in this case as = |cra| 
and (Jeff{e) = cFa + \cra\- It gives 

This means that a behaviour of ai,ff{e) in the percolation theory is completely 
determined by its odd part. From a general discussion of the functional struc- 
ture of the conductivity of the two-dimensional randomly inhomogeneous sys- 
tems and the known experimental and numerical results it is known that in 
the percolation limit the effective conductivity Uef / must have a nonanalytical 
behaviour near the percolation edge a;c = 1/2 or at small e > 

aeff{e) ^ (Ji{x - ^ ait\ (45) 

where a critical exponent of the conductivity t is slightly above 1 and can be 
represented in the form t — 1 + 5. Since the values of this exponent found by 
the numerical calculations are confined to be in the interval (1,10 - 1,4) [14], 
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then 6 have to be small and belongs to the interval (0.1 - 0.4). It follows from 
general formula (44) that one must have 

/a(e,^) (46) 

It means that the function <I>(e, z) has at small e some crossover on z under 
z ^ 1 from a regular (analytical) behaviour to a singular one. At the moment 
an exact form of this crossover is unknown. For example, it can be of the form 

<i>(e,z)^(<i>o(z2) + $i(zV)''''^', ^-1, (47) 

where $o(-z^) — > 0, $1(2:^) ^ $1 7^ and S{z) S ^ when z 1. 

However, as it follows from the formulas obtained above, one gets always 
(Jeff ^ in the limit 172 ^ , except the region near x = 1. It means 
that all these formulas obtained in FMSA approximation are not valid in the 
limit CT2 —> 0. This is confirmed by the appearance of the divergencies in the 
expansions of Te// in small concentrations in the limit z 1. This fact is 
a consequence of the made approximation. For example, in case of the model 
of the layered squares this is due to the "closing" or "locking" effect of the 
layered structure in the adopted approximation in the limit 172 —> 0. In order 
to obtain a finite conductivity in this model above threshold concentration Xc 
one needs to take into account the correlations between adjacent squares. It 
is easy to show that near the threshold an effective conductivity is determined 
by random conducting clusters formed out of the crossing random layers from 
neighboring elementary squares. As is well known, the mean size of these clusters 
diverges near the percolation threshold |l^ and for this reason the FMSA 
approximation cannot be applicable for the description of (Jeff hi the limit 
2: — » 1 (e > 0) and of the percolation problem. It follows from figs. 3 that the 
EM approximation overestimates Ue/ / , and both other formulas underestimate 
it in the region z — )■ 1, e > 0. We hope to investigate this limit in detail in the 
subsequent papers. 

Thus we have discussed the possible functional forms of the effective con- 
ductivity of the two- phase system at arbitrary values of concentrations. A new 
functional equation, generalizing the duality relation, was deduced in the FMSA 
approximation and its solution was found. We have constructed also a hierar- 
chical model of the random inhomogeneous medium and have found its effective 
conductivity in the same approximation at arbitrary phase concentrations. All 
formulas for the effective conductivity have the different functional forms. They 

(1) satisfy all symmetries including the dual symmetry and all necessary 
inequalities, 

(2) reproduce the general formulas for Ue// in the weakly inhomogeneous 
case. 

All these results confirm a conjecture that, in general, cTe// of the two- 
phase randomly inhomogeneous system may be a nonuniversal function and 
can depend on some details of the structure of the randomly inhomogeneous 
regions. Analogous conclusions were done in paper fl^ during a discussion of 
the possibility to find the generalization of the formula (3) for case > 3 . 
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The obtained formulas can be used for the approximation of the effective 
conductivity of some real randomly inhomogeneous systems like a corresponding 
formula of the EM approximation. The introduced composite method of the 
construction of the model random medium can be generalized on the other 
ways of determination of the effective intermediate conducting boxes. It can be 
done for different types of boxes as well as for different numbers of the possible 
types of the boxes. It is clear that then one will have to use another formulas 
instead of (3). 

This work was supported be RFBR grants 00-15-96579 and 02-02-16403. 
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